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CO ■ We report two new low metallicity blue compact dwarf galax- 

ies (BCDs), WISEP J080103.93+264053.9 (hereafter W0801+26) and 
WISEP J170233.53+180306.4 (hereafter W1702+18), discovered using the 
■ Wide-field Infrared Survey Explorer (WISE). We identified these two BCDs 

from their extremely red colors at mid-infrared wavelengths, and obtained 
Oh 1 follow-up optical spectroscopy using the Low Resolution Imaging Spectrom- 

O ■ eter on Keck I. The mid-infrared properties of these two sources are similar 

to the well studied, extremely low metallicity galaxy SBS 0335-052E . We de- 
termine metallicities of 12 + log(0/H) = 7.75 and 7.63 for W0801+26 and 
W1702+18, respectively, placing them amongst a very small group of very 
metal deficient galaxies (Z < 1/10 Z@). Their > 300 A H/3 equivalent widths, 
similar to SBS 0335-052E, imply the existence of young (< 5 Myr) star form- 
ing regions. We measure star formation rates of 2.6 and 10.9 M yr _1 for 
W0801+26 and W1702+18, respectively. These BCDs, showing recent star 

q . formation activity in extremely low metallicity environments, provide new 

laboratories for studying star formation in extreme conditions and are low- 
redshift analogs of the first generation of galaxies to form in the universe. 
Using the all-sky WISE survey, we discuss a new method to identify similar 
star forming, low metallicity BCDs. 
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Introduction 



In the local universe, we observe a class of galaxies known as blue compact dwarf 
galaxies (BCDs) which show signatures of strong, recent and ongoing star-formation activ- 
ity in low metallicity environments. These galaxies, which appear to be experiencing their 
first major burst of star formation, are characterized as generally having compact sizes 
(< 1 kpc), low luminosities (Mb > —18), strong emission line spectra, blue optical colors 
(U — B ~ —0.6), and their metal distribution peaks at one-tenth solar [Z = 1/10 Z Q ; 



Kunth fc 6stlinll2000l ). Typical H/3 equivalent width s range from 30 to 200 A, implying 
a burst of star formation no older than a few Myr kunthl[l999l ). BCDs offer pristine 
laboratories from which to study star formation processes, enrichment mechanisms for 
the interstellar medium, and the star formation history of galaxies in low metallicity and 
extremely dusty environments. In particular, BCDs are likely good analogs of the first 
generation of galaxies to form in th e universe, but are at redshifts accessible to detailed 
study (e.g., Meier fc Terlevichlll98lh. For a com prehensive review of BCDs and the most 
metal-poor galaxies, see iKunth fc Ostlinl ( 120001 ). 



Although discovered as early as the 1960's, BCDs and extremely low metallicity 
galaxies have remained an elusive population. The difficulty in finding low metallicity 
BCDs can be viewed as stemming from several factors. First, they have low luminosities, 
generally much fainter than L*. Second, at the lowest metallicities, below 0.1 solar, the 
major coolant species shift from [O III] to H and He ( IKunth fc Sargentlll986l ). resulting 
in lower [O III] A5007/H/3 ratios. This made s uch sources less readily identified in the 
previous generation of objective prism surveys (IKunth fc Sargentl Il986). though recent 
fiber spectroscopy surveys have less significant biases. Finally, they are rare. The most 
abundant samples of BCDs have been acquired by the Sloan Digital Sky Survey (SDSS) 
which has spectroscopic ally identified hundreds of BCDs at redshifts of a few tenths 



(e.g. 



Izotov et al.l 120061 ). implying a surface density of 0.1 BCD per deg 2 to the limits 
of the SDSS. However, SDSS rar ely finds the most ex treme metal poor systems, galaxies 
with Z < 1/10 Zq. For example. iKniazev et al.l ( 120031 ) estimate 0.004 deg -2 for extremely 
metal poor galaxies [XMPGs; 12 + log (O/H) < 7.65] for r < 17.77 mag. 



The two lowest metallicity BCDs known to date are SBS 
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20041 IWu et al.l 120071 ) . both of which have been extensively studied across the electro- 
magnetic spectrum. Though the metallicities for these two systems are very similar, 
12 + log (O/H) < 7.3 or Z < 1/40 Z Q , these galaxies exhibit very different morphological 
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and physical properties. iHirashita fc Huntl (120041 ) investigated these differences and sug- 
gest that they are best understood if two different modes of star formation are considered, 
an "active mode" and a "passive mod e". SBS 0335-052E is an active-mode BCD while 
I Zw 18 is a passive-mode BCD. In the IHirashita fc Huntl (120041 ) framework, active BCDs 
host super star clusters (SSCs) and are characterized by compact star formation regions 
(radius < 50 pc), high gas densities (> 500 cm -3 ), rich H 2 content, large dust optical 
depth, and high dust temperature. Passive BCDs, in contrast, have more diffuse star 
forming regions (> 100 pc), lower gas density (< 100 cm -3 ), cooler dust, and lack SSCs 
and large amounts of H 2 . 



Hirashita fc Huntl (J2004J) describe how the compactness of the star forming regions 



can simultaneously explain most of the observed differences between BCDs. Compact 
star-forming regions will rapidly become dusty from Type II supernovae, and this dust 
will be effective at reprocessing photons into the infrared. In particular, the gas free-fall 
time scale in such systems is less than 5 Myr, which leads to run-away star formation and 
the efficient creation of large quantities of dust. In contrast, the dynamical time of diffuse 
star forming regions, characteristic of passive BCDs, is longer than 10 7 yr. Such regions 



thus have lower star formation rates and are less infrared luminous. Hirashita fc Hunt 



(120041 ) suggest that the physical state of the gas might be driven by the size and spatial 
distribution of dust grains. 



Wu et ajj (120081 ) show that the infrared luminosities of these prototypical BCDs differ 
by two orders of magnitude, with SBS 0335-052E being ~ 100 times more infrared lumi- 
nous than I Zw 18. The mid-infrared emission of similar low metallicity BCDs originates 
mainly from hot (200 - 1500 K) dust, causin g active BCD s to ha ve significantly red colors 
at mid-infrared wavelengths. For exa mple, 



Houck 



of SBS 0335-052E peaks at ~ 28 /im. iDale et al.l (l2001bl ) model the infrared emission of 



et al. 



(120041 ) finds that the spectrum 



SBS 0335-052E and conclude that it can be characterized as having two dust components, 
a warm (~ 80 K) and a hot (~ 210 K) component. The prevalence of hot dust in active 
BCDs suggests that searching for the most extreme low metallicity, star-forming BCDs 
at mid-infrared wavelengths may be fruitful. 

This Letter reports the discovery of two n ew low metallicity BCDs discovered by 
the Wide-field Infrared Survey Explorer (WISE; IWright et al.ll2010l ). Launched on 2009 
December 14, WISE completed its first coverage of the entire sky on 2010 July 17, ob- 
taining at least 8 exposures per sky position in four passbands, 3.4, 4.6, 12 and 22 /j,m 
(W1,W2,W3 and W4). The corresponding 5cr point source sensitivities in unconfused 
regions are better than 0.08, 0.11, 1 and 6 mJy, respectively, and the point spread func- 
tion FWHM values for the four WISE bands are 6'.'1, 6'.'4, 6'.'5 and 12'.'0, respectively 
(IWright et al.ll2010l ). Focussing on two early examples of VF/S'i^-identified, low metallicity 
BCDs, we discuss the physical properties of these two new BCDs and place them within 
the context of recent BCD literature. Tsai et al. (in prep.) presents a larger sample of 
several dozen newly identified BCDs and discusses the selection technique in more detail. 
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We assume a solar metallicity (Z Q ) of 12 + log (0/H) = 8.96 flAllende Prieto et al.ll200lh 
in this paper. Magnitudes are reported in the Vega system and, where necessary, we 
adopt the concordance cosmology with Qm = 0.3, Q\ = 0.7 and Hq = 70 kins -1 Mpc -1 . 



2. BCD Selection and Keck Spectroscopy 

We identified WISEP J080103.93+264053.9 (hereafter W0801+26) and WISEP J170233.53+180. r 
(hereafter W1702+18) from their extreme mid-infrared colors and relatively bright optical 
flux densities (B < 19.5, from the USNO-B1 catalog). The mid-infrared photometry was 
acquired from the WISE Preliminary Release Source Catalog (thus the "P" designation 
in the WISE nomenclature). The simultaneous multi-wavelength WISE observations of 
W0801+26 and W1702+18 were carried out on UT 2010 April 12 and UT 2010 March 2, 
respectively, and provided total exposure times of ~ 100 and 120 sec, respectively. For 
details on the WISE processing in the Preliminary Release, see the WISE Preliminary 
Release Explanatory Supplement^. WISE images of the two new BCDs are presented in 
Fig. [TJ For comparison, we also show WISE images of the two prototypical low metall- 
city BCDs, SBS 0335-052E and I Zw 18. Photometry for all four sources is presented in 
Table 1. The new BCDs have unique and easily identified mid-infrared colors, similar to 
SBS 0335-052E . In contrast, I Zw 18 is unremarkable in the WISE data. None of the 
sources are resolved by the 6" WISE beam. The signal-to-noise ratios of SBS 0335-052E 
and the new BCDs are > 25 in all four WISE bands. 

Noting their unusual WISE colors and relatively bright optical magnitudes, we ob- 
tained optical spectroscopic follow-up observations of W0801+26 and W1702 +18 during 



twilig ht using the dual-beam Low Resolution Imaging Spectrometer (LRIS; lOke et al. 



19951 ) on the Keck I telescope. We observed W1702+18 on UT 2010 March 12 and 
W0801+26 on UT 2010 November 8. The conditions were photometric on both nights, 
and both observations used the l'/5 wide longslit, the 5600 A dichroic and the 400 
£mm _1 grating on the red arm of the spectrograph (blazed at 8500 A; resolving power 
R = A/AA ~ 700 for objects filling the slit). The March 2010 observations used the 
600 £mm _1 grism on the blue arm (blazed at 4000 A; R ~ 1000), while the November 
2010 observations used the slightly lower resolution 400 £mm~ 1 grism (blazed at 3400 A; 
R ~ 600). A single 300 second exposure was obtained of each source during the observing 
runs. We processed the data using standard procedures , flux calibrated the spectra using 



observations of standard stars fro mlMassev Gronwal 



extinction using the dust maps of ISchlegel et al.l ( 119981 ). 



(Il990l ). and corrected for Galactic 



The final, reduced spectra show a multitude of extremely high equivalent width, 
narrow emission lines at low redshift and are presented in Fig. |2j We have used the ratio 
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of Ha to H/3 and H7 to H<5 emission line strengths to derive and cor rect for intrins i c dust 
extinction assuming Case B recombination and the extinction law of ICardelli et al.l (Il989f ) 
with Ry = 3.1. For W0801+26, the Balmer line ratios are close to the extinction-free 
ratios of Case B, implying no significant extinction is observed. For W1702+18, we find 
an intrinsic extinction of Ay = 0.18. The derived redshifts are presented in Table 1, while 
Table 2 presents the extinction-corrected line flux measurements and errors obtained using 
the SPECFIT package within IRAF. Equivalent widths from H/3 are also given in Table 2. 



3. Results 
3.1. Mid-Infrared Properties 

Fig. [1] shows WISE color images of W0801+26 and W170 2+18, as well as of proto 



typical BCDs SBS 0335-052E and I Zw 18. In the framework of iHirashita fc Huntl (I2004T ). 
SBS 0335-052E is an active BCD, while I Zw 18 is passive BCD. Most Galactic stars and 
low redshift galaxies are dominated by a Rayleigh- Jeans spectrum at the wavelengths ob- 
served by WISE, providing blue colors in the color mapping shown in Fig. [TJ In contrast, 
active BCDs contain substantial quantities of dust in a dense cloud, which give them 
unusual and readily-identified red colors at near- and mid-infrared wavelengths. BCDs 
like I Zw 18, on the other hand, do not contain hot dust, and thus are undistinguished in 
the WISE bands (see Table 1). 

As discussed in Yan et al. (in prep.), sources with red colors across the first two WISE 
bands, Wl — W2 > 2, are relatively rare. For example, the COSMOS field has only one 
source with such red WISE colors in a field slightly larger than 1 deg 2 (Stern et al., in 
prep.). The coolest brown dwarfs, with spectral types later than T5, have red Wl — W2 



color s due to methane abso rption at 3.3 /im (e.g., iBurgasser et al.ll2011uKirkpatrick et al. 



20 lit iMainzer et al.l 120111 ) and are quite faint at optical wavelengths. The WISE team 
has also been pursuing a population of galaxies with red W2 — WA colors which appear to 
be dominated by heavily obscured AGN at z ~ 2 (e.g., Eisenhardt et al., in prep.). Such 
galaxies are also quite faint optically. In contrast, as shown in Table 1, low redshift BCDs 
have extremely red WISE colors but are bright at optical wavelengths, providing for their 
easy identification and follow-up. Tsai et al. (in prep.) discusses the WISE selection of 
BCDs in more detail. 



3.2. Spectroscopic Properties and Metallicity 

The Keck spectra of W0801+26 and W1702+18, presented in Fig. El show a multitude 
of strong hydrogen and oxygen emission lines as well as weak stellar continuum. These 
spectra are characteristic of sources whose optical emission is dominated by H II regions. 
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Key line measurements for W0801+26 and W1702+18 are presented in Table 2. Standard 
diagnostic diagrams, such as comparing [0 III]/H/5 to [N II] /Ha, imply that the emission 
lines for both new Wiff-E-se lected BCDs are driven by star formati on processes rather 
than by AGN activity (e.g., baldwin et allll98ll : iKewlev et aDboOlh . 



One of the most striking features of the spectra presented in Fig. [2] is the weakness 
of the [N II] lines relative to the Ha line that they flank. This immediately s uggests that 



both g alaxies have extremely low metallicities. Following the prescriptions of llzotov et al. 
(booq ). we derive electron temperatures T e ([0 III]) ~ 1.7 x 10 4 K for W0801+26 and 



- 1.9 x 10 4 K for W1702+18, from their [O III] line strengths at 4363, 4959 and 5007 
A. Based on the relative strengths of the Balmer lines to these oxygen lines, we then 
derive total heavy element abundances and metallicities of 12 + log (O/H) = 7.75l^i7 an d 
12 + log (O/H) = 7.63 ± 0.06 for W0801+26 and W1702+18, respectively, where these 
errors were derived using a Monte Carlo method with 10,000 simulations. 

Using th e Ha and 22 ^m flu x densities we measure star formation rates (SFRs) as 
prescribed in lCalzetti et~al1 fl2007h . We find SFR = 2.6 and 10.9 M yr" 1 for W0801+26 
and W1702+18, respectively, which imply that these two BCDs are currently undergoing 
very extreme bursts of star formation activity. In addition, the H/3 equivalent widths > 
300 A for W080 1+26 and W1702+18 im ply that the ages of the stellar populations are 
less than 5 Myr (ISchaerer fc Vaccalll998l ). 



4. Discussion and Summary 



The low metallicity and mid-infrared dominated spectral energy distributions of these 
newly identified, WISE-selected galaxies place them in a rare class of active, low metallic- 
ity BCDs. We compare the infrared properties of these BCDs with the infrared properties 
of known BCDs and find that the majority of low metallicity BCDs in the literature are 
several magnitudes fainter at 12 fim (W3) than the two new sources. The properties of 
the two new BCDs are strikingly similar to one of the most well known BCD, SBS 0335- 
052E, suggesting that the mechanisms responsible for SBS 0335-052E might be the same 
mechanism at work in W0801+26 and W1702+18. 

In Figure 3 we show WISE Wl — W2 color versus metallicity, 12 + log (O/H), for 
different samples of BCDs. After v isual inspections to remove g alaxies that c l early are not 
BCD s, we include U BCD s from IWu et all (hoosh . 173 from llzotov et al.l tood) thre e 
from lKniazev et all (l2003h. SBS 0335-0 52E, I Zw 18, HS 0822+3542 (llzotov et al.ll2006h . 
HS 2236+1344 rtlzotov fc Thuanl 12007^ and our two newly discovered BCDs. With the 
exception of I Zw 18 and HS 0822+3542, we observe a slight correlation between W1 — W2 
and metallicity for these samples, with Wl — W2 color becoming increasingly redder with 
decreasing metallic ity. We apply the Spearman's rank correlation coefficient test on the 
Izotov et al.l (120061 ) sample and recover a Spearman coefficient of p = —0.36. For a sample 
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of that size, this indicates a noticeable correlation, deviating from the null-hypothesis at 
the 4<T level. 

Since polycyclic aro matic hydrocarbon (PAH ) emission is largely absent in low- 
metallicity environments (jEngelbracht et al.l 120081 ) . the mid- infrared spectra of BCDs 
i s dominated by therma l emission with a minor contribution from nebula continuum 
dSmith fc Hancock! hood ). Thus, the 3.4 /im to 4.6 /zm flux ratio is a relatively clean tem- 
perature indicator for the hot dust heated by young stars. At the lowest metallicitities, 
12 + log (O/H) < 7.8, BCDs seem to group into two different color regions, W1 — W2 ~ 2.0 
and ~ 0.5. This suggests two distinct star-formation processes can be d o minan t, support- 
ing the passive and active BCD modes suggested by 



Hirashita fc Huntl (120041 ) . 



We are currently conducting follow-up observations of many more promising can- 
didates discovered using the extreme infrared properties exhibited by these two sources 
and SBS 0335-052E. These results will be presented in a follow-up paper (Tsai et al., in 
prep.). We expect that WISE will significantly increase the number of low metallicity 
BCDs known, perhaps finding the most pristine example in the local universe. Such stud- 
ies will allow us to construct a statistical sample from which to study and characterize the 
BCD population. The WISE BCDs will be a unique sample of mid-infrared bright, low 
metallicity galaxies for studying dust grain formation in low metallicity environments, and 
will help us to understand the thermal dust emission of high-redshift, starburst galaxies. 
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Fig. 1. — WISE color images (blue = Wl, green =W2, red = W3) for the two newly 
discovered low metallicity BCDs as well as the prototypical BCDs SBS 0335-052E and 
I Zw 18. Images are 5' on a side, with North up and East to the left. "Active" BCDs 
such as SBS 0335-052E and the two new BCDs have unique, extremely red WISE colors, 
making them readily identifiable from the WISE all-sky survey. 
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Fig. 2.— Spectra of W0801+26 and W1702+18, obtained with the LRIS instrument 
on Keck I. Both spectra show strong, narrow, high equivalent width emission lines from 
hydrogen and oxygen, as well as weak [N II] flanking the Ha emission line (insets). The 
line ratios imply that the emission is powered by star formation, not AGN activity, from 
very low metallicity systems. 
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Fig. 3. — Mid-infrared color (Vega) vs. metallicity, 12 + log (O/H), for several samples of 
BCDs, as indicated. The plotted sources are required to have SNR > 5 in Wl and W2. For 
comparison, we also plot SBS 0335-052E, I Zw 18, HS 0822+3542 and HS 2236+1344, 
four of the lowest metallicity galaxies known to date. W0801+26 and W1702+18 are 
significantly redder in Wl — W2 than the majority of known BCDs. 
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Table 1. Properties of Low Metallicity BCDs. 





W0801+26 


W1702+18 


nrjn noocr ntro"!'^ 

bnb 0335-052L 


T 7 -1 o 

1 Zw 18 


R.A. (J2000) 


08:01:03.929 


17:02:33.534 


03:37:44.032 


09:34:02.088 


Dec. (J2000) 


+26:40:53.91 


+ 18:03:06.44 


-05:02:40.26 


+55:14:26.23 


z 


0.0260 


0.0425 


0.0134 


0.0025 


M B (mag) 


-16.8 


-18.3 


-17.1 


-15.4 


B (mag) 


18.18 


18.39 


16.48 


14.46 


Wl (mag) 


15.03 ± 0.04 


14.26 ± 0.03 


14.53 ± 0.03 


15.39 ± 0.04 


W2 (mag) 


12.94 ± 0.03 


12.03 ± 0.02 


12.52 ± 0.03 


14.86 ± 0.07 


W3 (mag) 


8.29 ± 0.03 


7.58 ± 0.02 


7.66 ± 0.02 


11.88 ± 0.14 


WA (mag) 


5.38 ± 0.03 


4.96 ± 0.03 


5.04 ± 0.03 


7.84 ± 0.14 


12 + log (0/H) 


7 7 cr+0.04 
' • ' d -0.17 


7.63 ± 0.06 


7.11-7.3 


7.17 



Table 2. Spectroscopic Line Measurements. 





W0801+26 


W1702+18 


[0 II] 3727+3729 


234.1 ± 2.6 


334.2 ± 11.1 


R5 4101 


69.5 ± 1.1 


127.9 ± 2.0 


H 7 4340 


132.5 ± 4.1 


232.1 ± 13.7 


[0 III] 4363 


41.5 ± 3.1 


84.2 ± 8.7 


H/3 4861 


247.3 ± 3.0 


438.5 ± 5.7 


[0 III] 4959 


527.8 ± 11.3 


921.8 ± 31.6 


[0 III] 5007 


1596.4 ± 18.9 


2679.2 ± 62.3 


Ha 6563 


657.4 ± 9.0 


1227.2 ± 27.5 


[Nil] 6584 


10.0 ± 0.8 


21.1 ± 0.6 


-EW(H/3) 


370.2 ± 6.1 


324.4 ± 6.4 



Note. - - The top rows present the line fluxes in 
units of 10~ 16 erg cm -2 s _1 . The final row gives the 
H/3 equivalent width (EW) in units of A. 



